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INTRODUCTION 


The  focus  of  this  research  is  to  develop  predictive 
techniques  for  determining  the  surface  temperature  of  bridge-type 
structures  using  known  or  assumed  weather  data*  This  technique 
would  be  used  in  Tactical  Decision  Aide  (TDA)  models.  TDAs  are 
used  to  provide  military  field  commanders  information  about  the 
apparent  contrast  between  a  target  and  its  background. (Goldman  et  al 
1983,  p.4^.  Having  this  information  before  deployment  increases 
the  effectiveness  of  infrared  guided  missiles,  or  helps  the  field 
commander  make  an  informed  decision  about  the  type  of  weapons  to 
use. 

As  used  in  military  research  or  operations,  a  TDA  consists  of 
three  major  segments  (Ibid).  These  are  prediction  of  inherent 
target-to-background  contrast,  atmospheric  effects  on  the  inherent 
contrast,  and  sensor  response  to  Che  apparent  contrast.  The  final 
result  is  a  predicted  acquisition  and  lock-on  range  for  a 
specified  target,  background,  and  weather  history.  This  research 
concentrates  on  the  first  of  these  segments. 

Contrast  predictions  of  this  type  have  been  successfully 
performed  by  personnel  of  the  Air  Force  Armament  Laboratory 
(AFATL)  (Goldman  et  al,  1983).  *^The  model  for  this  effort  was  a 
one-dimensional  Runge-Kutta  analysis  of  the  heat  flow  equations. 
However,  these  predictions  were  typically  for  small  objects  such 


as  vehicles  or  concrete  bunkers*  Modelling  large  objects  such  as 
bridges  presents  certain  difficulties  as  noted  by  Goldman  and 
coworkers  ^(1983,  pp • 34-38 ) , 

Another  related  work  of  this  type  was  by  Thepchatri  and 
coworkers  in  1977*^  This  effort  concentrated  on  predicting  the 
internal  temperature  profile  of  a  bridge  for  the  purpose  of  stress 
analysis.  However,  the  method  used  is  directly  related  to  this 
effort  because  the  Internal  temperature  profile  is  dependent  on 
the  surface  temperature. 

Thepchatri  and  coworkers  used  two  types  of  models;  a 
one-dimensional  finite  difference  form  of  the  heat  flow  equation, 
and  a  two  dimensional  finite  element  analysis  of  the  heat  flow. 


I.  APPROACH  TO  TH£  PHYi>ICAL  PKOBL£H 


The  basic  efforc  here  was  to  deteralne  specltlc  factors 
Involved  In  predicting  surface  temperatures  of  bridges  and  develop 
an  adequate  approach  for  doing  so.  AFATL  had  taken  radiance 
oeasureiaents  of  a  bridge  In  the  Oestln,  KL  area,  so  it  was  decided 
to  use  this  structure  to  evaluate  modelling  techniques. 

The  Oescln  bridge,  shown  In  Figure  I,  runs  east -west  over  a 
bay.  This  provides  a  constant  background  which  Is  simpler  to 
model.  Also,  because  It  is  an  east-vest  bridge,  heating  by  solar 
radiation  takes  place  almost  exclusively  on  the  top  surface 
(Thepchatrl  et  al,  1977,  p.40).  Since  solar  radiation  is  the 
dominant  energy  exchange  mechanism,  a  one-dlmenslonal  model  was 
thought  to  be  adequate  to  predict  Che  roadway  surface  temperature. 

This  hypothesis  Is  further  substant iated  by  thermograph 
charts  of  the  Destln  bridge  (Appendix  A)  which  show  the  upper 
surface  of  thr  roadway  as  the  dominant  feature,  exhibiting  the 
greatest  tt^mperature  change  throughout  a  day's  cycle.  These 
thermograph  charts  are  from  measurements  made  by  AFATL  personnel 
on  17  May  1984  (TABILS  1).  The  measurements  were  made  with  a  AGA 
680  Thermovision  camera  in  two  wavelength  bands:  one  in  the  3.0 
to  5.0  micrometer  region  and  the  other  In  the  7.0  to  14.0 
micrometer  region.  The  charts  in  Appendix  A  are  scenes  in  the 


approxlaately  the  eceae  by  (her»uvl:>lon  cattera. 


Goldaan  aad  co^ork^era  (t9d3,  p.33)  noted  chac  tot  atroctufc^a 
such  as  bridges,  cargec  recognlrloa  ney  be  accaapHshed  by  virtue 
of  coQCrasts  wlchla  the  target.  Therefore,  In  order  to  study  the 
contrasts  within  the  bridge,  it  was  thought  necessary  to  provide 
soae  Mans  of  predicting  the  t  «perature  of  these  parts  of  the 
bridge  in  addition  to  the  roadway  structure. 


The  thertsograph  charts  indicate  that  the  sides  of  the  roadway 
and  the  supporting  structure  teaperatures  tend  to  follow  each 
other.  For  simplicity*  it  was  decided  to  model  Just  the  vertical 
piling  support  underneath  the  roadway  structure  and  assume  that 
this  temperature  represented  the  temperature  of  the  sides  also.  A 
schematic  of  the  bridge  as  modelled  is  shown  in  Figure  2. 


Modelled  in  this  fashion,  the  piling  model  is  independent  of 
the  model  for  the  roadway  surface,  which  neglects  conduction  and 
radiation  between  the  piling  and  the  roadway  structure.  This  is, 
however,  consistent  with  the  conclusions  of  Goldman  and  coworkers 
(1983,  p.36)  that  bridges  can  be  separated  into  roadway  and 
supporting  structures. 
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A  one-dimensional  model  similar  to  Chat  used  by  AFATL  Co 
model  vehicles  was  used  Co  model  Che  DesCin  bridge,  and  Chese 
resulcs  were  compared  Co  cemperaCure  daca  derived  from  Che  17  May 
1984  measuremencs*  However,  noC  all  bridges  will  exhibic  chls 
simple  geomecry,  and  for  bridges  running  a  direccion  ocher  chan 
easc-wesc,  solar  radiacion  on  Che  sides  should  be  accounCed  for 
using  a  Cwo-dimenslonal  model  (ThepchaCri  ec  al,  1977,  p.40).  For 
chls  reason,  a  cwo-dimensional  model  was  used  firsc  Co  model  che 
Descln  bridge  in  ics  accual  conf iguracion  for  comparison  wlch 
resulcs  of  Che  one-dimensional  model.  Then  che  Descln  bridge  was 
modelled  as  if  ic  were  ac  a  norch-souch  oriencacion  Co  decermlne 
Che  effecc  of  solar  loading  on  che  sides. 
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II.  MODEL  DEVELOPMENT 


The  models  used  in  earlier  studies  discussed  previously,  were 
based  on  the  Fourier  heat  conduction  equation  (Kreith,  1973,  p.84). 


3T  9^T  3^1  3^1 

-  -  K(  --  t  --  +  -5  ) 
3t  9x^  3y^  3z 


where  t  refers  to  time;  x,  y,  and  z  refers  to  the  spatial 
Cartesian  coordinates;  and  K  is  the  thermal  diffusivity.  This 
assumes  there  are  no  internal  heat  sources  and  that  the  material 
is  isotropic.  The  equation  is  subjected  to  the  boundary  condition 
of  a  heat  flux  at  the  surface  dependent  on  the  environmental 
conditions. 


Environmental  Interactions: 


As  mentioned  above,  solar  radiation  is  the  dominant  energy 
exchange  mechanism  affecting  surface  temperature.  However, 
surface  temperature  is  also  Influenced  by  ambient  radiation, 
convection,  and  conduction  to  the  interior  surfaces  (Goldman  et  al 
1983,  pp.5-6  and  Thepchatrl  et  al,  1977,  p.23),  with  conduction 
being  accounted  for  by  the  heat  conduction  equation. 


heat  transfer  mechanism  for  concrete  is  zero  (Goldman  et  al,  1983, 
p«12).  Also,  evaporation  and  condensation  are  not  considered  in 
the  models  developed  by  Thepchatri  and  coworkers.  For  these 
reasons,  evaporation  and  condensation  were  not  considered  in  this 
effort. 

Solar  and  ambient  radiation  and  convection  on  a  bridge 
surface  are  dependent  on  the  local  weather  conditions.  Actual  use 
of  a  model  such  as  this  in  an  operational  military  environment 
will  likely  have  to  rely  on  predicted  or  estimated  weather 
conditions.  Both  Goldman  and  Thepchatri  and  their  coworkers 
(Goldman  et  al,  1983,  pp. 19-27  and  Thepchatri  et  al,  1977, 
pp. 14-18,  23-26)  discuss  methods  of  making  predictions  of  the 
radiation  quantities. 

For  this  study,  weather  data  taken  in  conjuction  with  the 
AFATL  measurements  of  the  Destin  bridge  was  used.  This  data  was 
taken  over  a  time  period  approximately  from  six  o^clock  in  the 
morning  until  midnight  on  the  day  of  the  measurements.  The  data 
was  taken  at  approximately  one  minute  intervals  with  the  following 
quantities  being  measured:  wind  speed,  wind  direction,  air 
temperature,  relative  humidity,  barometric  pressure,  solar  and 


ambient  radiation. 


The  weaCher  daca  Indicates  chat  the  wind  was  froa  a  generally 
northerly  direction  and  the  velocity  was  noraally  10  knots  or 
less,  but  occasionally  as  high  as  I?  knots.  The  relative  humidity 
varied  froa  about  20X  to  as  high  as  45t,  and  the  barometric 
pressure  varied  froa  about  1021  to  1024  millibars.  The  variation 
in  air  temperature  throughout  the  day  Is  shown  In  Figure  3,  and 
Che  variation  of  Che  solar  and  ambient  radiation  Is  shown  In 
Figure  4. 

The  solar  and  ambient  radiation  quantities  were  measured  by 
means  of  a  pyrhellometer  and  a  pyranometer  respectively.  A 
pyranomeCer  Is  a  device  which  measures  the  total  radiation  within 
Its  hemispherical  field  of  view  with  the  direct  beam  radiation 
blocked  from  the  sensor  surface  by  an  occulting  disk  (Krelch  and 
Krelder,  1978,  pp. 63-66).  A  pyrhellometer  Is  a  similar  device 
which  measures  the  direct  beam  solar  radiation  on  a  horltzontal 
surface  while  blocking  the  diffuse  sky  radiation  by  mounting  the 
detector  at  the  base  of  a  Cube  pointed  directly  at  Che  sun  during 
Che  day  (Ibid). 

When  the  sun  Is  shining  on  the  bridge  surface,  the  heat 
absorbed  by  the  horltzontal  surface  Is  given  by  Che  Intensity  of 
solar  radiation,  as  measured  by  the  pyrhellometer, 

multiplied  by  the  solar  absorptivity  of  the  surface.  The  equation 
(Thepchatri  et  al,  1977,  p.23)  is 


The  solar  radlaclon  absorbad  Is,  of  course,  ac  a  very  shore 


wavelength.  Radiated  heat  exchange  at  longer  wavelengths  Involves 
a  balance  between  Che  aeblent  radiation,  U^ky.  measured  by  the 
pyiTduometer  and  the  eaXtted  radiation  ot  the  surface  which  follows 
the  Stefan-Boltzaann  Law  (Ibid,  p.24).  So,  given  an  average 
emlsslvlty,  e,  of  the  surface  the  heat  loss  at  the  top  surface  'o 
the  sky  Is  given  by 


Ql  •  elo(T^^273.l5)^-U^^y) 

where  Tg  Is  the  surface  teaperature  of  the  bridge  In  degrees 
Celsius  and  0  is  the  Stefan-Boltzaann  Constant. 

However,  the  pyranoaeter  aeasureaent  will  Include  the 
contribution  of  clouds,  which  absorb  radiation  during  the  day  and 
emit  It  back  to  the  earth  during  the  night  (Ibid,  p.26). 

Therefore,  for  surfaces  which  do  not  face  the  sky,  an  estimation 
of  the  ambient  radiation  which  does  not  Include  the  cloud 
contribution  Is  required.  This  Is  estimated  by  applying  the 
Stefan-Boltzmann  Law  to  the  ambient  air  temperature,  T^.  Thus, 
the  heat  loss  to  the  surroundings  from  surfaces  not  facing  the  sky 
Is  given  by 


Ql  -  ea[(Tg  +  273.15)^  -  (T^  +273.i:)''j. 
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Another  major  environmental  interaction  is  convection, 
including  both  free  and  forced  convection.  At  the  bridge  surface, 


this  is  given  by 


where  is  the  wind  speed  and  hj^  and  h2  are  empirically  derived 
constants  (Ibid,  p.27)  which  are 


hj  -  3.775 

h2  -  1.689  W-s/m^-°C 


Another  interaction  of  the  bridge  with  the  environment  which 
sometimes  needs  to  be  considered  is  the  reflected  and  emitted 
radiation  of  the  water  below.  This  is  a  factor  which  Goldman  and 
coworkers  did  not  consider  because  they  did  not  consider  objects 
with  the  lower  surface  open  to  the  environment.  Thepchatrl  and 
coworkers  (1977,  p.61)  did  make  a  rough  estimate  of  the  reflected 
radiation  by  assuming  the  bottom  surface  absorbed' 10%  of  the 
incoming  solar  radiation. 


For  this  analysis,  it  is  assumed  that  the  water  reflects  the 
incoming  solar  radiation  and  the  ambient  radiation  and  that  this 
is  absorbed  by  the  bridge  surface  like  the  ambient  radiation  when 
appropriate.  This  is  given  by 


W  -  '  r- 


where  r  is  the  reflectivity  of  the  water.  Also  note  that  the 
average  emlsslvlty  used  for  the  Interaction  with  ch«  aablent 
radiation  is  used  here  as  the  absorptivity. 

The  water  also  emits  radiation,  and  this  is  also  absorbed 
like  the  ambient  radiation  given  as  follows 

Qg  -  ee„(o(T^  +  273.15)^] 

where  is  the  emissivlty  of  the  water  and  is  the  temperature 
of  the  water  in  degrees  Celsius. 

One-Dimensional  Model; 

The  simplest  approach  to  predicting  the  bridge  surface 
temperature  is  a  numerical  approximation  of  the  Fourier  equation 
in  one  dimension.  The  two  studies  which  have  been  discussed  here 
each  used  slightly  different  approaches  for  this  approximation. 
These  are  briefly  discussed  in  this  subsection.  However,  for  the 
purposes  of  this  study,  a  one-dimensional  model  was  developed 
which  combines  features  of  the  two  previously  developed  models 
best  suited  for  this  purpose  or  for  ease  of  application. 


The  model  developed  by  Goldman  and  coworkers,  as  previously 
mentioned,  is  a  Runge-Kutta  approximation  of  the  Fourier 
equation.  This  is  accomplished  by  first  approximating  the 
second  order  spatial  derivative  with  a  finite-difference 
technique.  This  yields  a  system  of  first  order  differential 
equations  in  time  which  are  solved  using  the  Runge-Kutta  technique 
(Goldman  et  al,  1983,  pp.6-11) 

The  system  of  equations  results  from  dividing  the  structure 
into  layers  for  application  of  the  finite-difference  technique, 
yielding  an  equation  for  each  layer.  In  this  case,  the  solution  is 
for  the  temperatures  at  the  centers  of  the  layers.  The  surface 
temperature  is  found  using  a  complicated  substitution  scheme 
outside  the  Runge-Kutta  loop  based  on  the  surface  temperature  of 
the  previous  time,  the  boundary  conditions,  and  the  first  layer 
temperature. 

The  study  by  Goldman  and  coworkers  did  not  discuss  modelling 
the  lower  boundary  open  to  the  air  as  is  the  case  with  a  bridge. 

This  can  be  easily  derived,  as  shown  later,  by  applying  the  same 
type  of  boundary  condition  on  the  last  layer  as  is  applied  on  the 
first. 

Thepchatri  and  coworkers  developed  a  full  finite-difference 
model  in  both  the  time  and  space  coordinates  (Thepchatri  et  al,  1977, 
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pp. 31-39).  In  this  approach,  the  structure  is  also  divided  into 
layers  but  with  the  temperature  solution  being  at  the  interface  of 
the  layers.  The  surface  temperatures  are  solved  as  part  of  the 
finite-difference  scheme  by  combining  the  equations  for  the  heat 
transfer  at  the  boundaries  with  the  finite-difference  equations  of 
the  top  and  bottom  layer  boundaries. 

The  model  developed  for  this  study  used  the  layer-centered 
temperature  solution  and  solved  the  surface  tempertures  as  part  of 
the  approximation  routine.  Figure  5  is  a  schematic  of  how  this 
type  of  structure  is  applied. 
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Figure  5.  Spacing  of  Temperature  Values  in  1-D  Model. 
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In  solving  for  the  temperature  at  the  center  of  a  layer  of 
thickness  d,  the  second-order  derivative  Is  approximated  using  a 
centered-dlf ference  form  as 


2 

ax  d^ 


Using  this  centered-dlf ference  form  keeps  the  calculation  of 
the  layer  temperature  within  the  layer.  Tliu!..  adjacent  laycr»  of 
different  materials  can  be  considered  without  calculating  average 
dlffuslvlty  values.  However,  except  for  the  surface  temperatures 
only  the  center  layer  temperatures  are  assummed  to  be  known. 

Thus,  the  temperatures  at  the  half-steps  are  assumed  to  be  the 
average  of  the  temperatures  of  the  two  adjacent  layers.  Applying 
this  Co  Che  Fourier  equaclon  gives 


- 2T,  .  Ti.i). 

dt  d2 


In  Che  Cop  layer,  the  temperature  at  is  T^.  So,  the 

approximation  to  the  Fourier  equaclon  at  Is 


dT  K 

--  -  -:;(T2  -  3T.  +  2T  ). 
dt  d^  ^ 


Slnllarly,  for  the  bottom  Nth  layer  where  is  Tl,  the 

approxlmatloa  Is 


dTj,  K 

—  -  -^(2Tl  -  3Tj.  + 
dt 


Applying  this  technique  to  the  top  surface  would  give 


dig  2K 

-  2T3  > 


But,  not  exist.  To  remedy  this,  the  heat  balance  at 

the  surface  Is  considered*  At  the  top  surface  of  a  roadway 
structure,  this  can  be  expressed  as  (Thepchatri  et  al,  1977,  p*33) 


Heat  absorbed  from  *  Heat  lest  by  convection  + 

short  wave  radiation  Heat  lost  by  long  wave  radiation  + 

Heat  lost  by  conduction. 


Algebraically,  this  is 


aWsun-  (hj  +  -  T^)  +  e(o(Ts  +  273.15)  -  Msky^ 


--<’■1  -  Vw) 


where  the  last  term  represents  the  conduction  at  the  surface 
expressed  in  centered-dif ference  form,  with  k  being  the  thermal 
conductivity.  Solving  this  expression  for  Tg_^j  and  substituting 
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that  into  the  approximatioa  of  the  Fourier  equation  at  the  top 
surface  gives 


dTg  AOc 

- -  "r(Ti  -  T2)  +  2W(Tg) 

dt  d^ 


where 


Wdg) 


--Wsun  +  --[Wsky  *  +  273. IS)"*] 

d  d 

+  (-^  +  -^Va)(T^  -  Tg). 
d  d 


The  alpha  terms  are  eavlronmental  constants  the  same  as 
derived  by  Goldman  and  coworkers  (1983,  pp.7-8).  These 


d  Pcd 

02  e 

d  Pcd 


^  .  -lli 

d  Ocd 


%  .  \ 

d  Dcd 


those 

are 


A  similar  scheme  is  used  at  the  lower  surface  of  the  roadway 


Heat  gained  by  conduction  Heat  lost  by  convection  + 

+  Heat  gained  from  water  «  Heat  lost  by  long  wave 
radiosity  radiation. 

Algebraically,  with  the  conduction  term  expressed  in  the  same 
manner  as  for  the  top  surface,  this  is 
k 

--'fu  -  ■'L4Jsd>  “sky)  +  +  273.15)'' 

-  '>>1+  -  ■'a) 

+  eoKTg  +  273.15)^  -  (I4  +  273. IS)*]. 

Solving  for  TL+ljd  substituting  into  the  Fourier  equation  gives 


-“(T„  -  T  )  +  2W(T  ) 
jj2  N  L  L 


where 


W(T,)  -  ~((Ta  +  273.15)^  -  (Tg  +  273.i5)^Ja  +  E 
d 

+  +  -\){Ta  -  T,)  7-  7-  ). 

d  d  d 


Additional  environmental  constants  in  this  case  are  as  follows: 


E„  -  ---(T  +  273.15)^0 
pcd 

“6 


It  is  important  to  note  here  that  the  emitted  and  reflected 


radiation  of  the  water  is  applied  only  to  the  lower  surface  as 
would  be  the  case  for  the  roadway.  Specific  applications  will 
determine  exactly  when  the  various  environmental  constants  should 
be  applied. 

This  system  of  expressions  for  the  Fourier  equation  can  be 
used  with  the  Runge-Kutta  method  to  obtain  solutions  of  the  time 
derivatives.  This  is  done  for  this  study  using  the  Runge-Kutta 
method  of  order  four  (Burden  and  Falres,  1985,  pp. 263-265).  In 
this  case,  the  equations  are 


An  initial  approximation,  T^q,  for  each  temperture  is 
assumed.  Then,  if  the  time  Increment  is  h  and  with  Tj^j 
representing  the  ith  temperature  at  the  Jth  time  Increment,  the 


Ich  temperature  at  the  next  time  Increment  is  obtained  by  first 
calculating 


hF^Ctj  ,T^j  , . . .  ) 


hFi(tj+-,Tsj+i5kis,  •  •  • 


hF i(t  j+^,Tgj+%k2s,  • . . 

hFi(t j+hjTgj+ksg, . . . »TLj+k3L). 


Then  the  solution  at  the  next  time  Increment  Is 


'^l,j+l  “  "^Ij  *  ^^11  2^21  ^*^31 


The  time  derivatives  can  also  be  expressed  In 
finite-difference  form,  which  gives  the  following  system  of 
equations : 


t;  -  Tg  +  4h-5(Ti-  Tg)  +  2hW(Tg) 
d 


T\  -  Tj  +  h-|(T2  -  3Ti  +  2Tg) 

d 


n  -  Ti  +  h-^(Ti+i  -  2Ti  +  Ti-i) 
d 


Tn  -  Tn  +  H<2Tl  -  3Tn  +  Tn_i) 


\  Tl>  *  ^hWCT^) 

d 

where  h  Is  again  the  finite  time  Increment  and  the  prime 
represents  the  new  temperature  at  the  end  of  the  time  step*  Both 
the  Runge-Kutta  Method  and  the  full  finite-difference  method  were 
used  In  this  analysis  for  comparison. 

Two-Dimensional  Model; 

A  two-dimensional  model  was  developed  to  study  what  effect 
solar  loading  on  more  than  one  side  might  have  on  the  surface 
temperature.  This  situation  is  significant  when  the  bridge 
Is  at  an  orientation  other  than  east-west  and  would  be 
particularly  significant  on  smaller  parts  of  the  bridge  such  as 
the  vertical  pilings. 

Thepchatrl  and  coworkers  (1977,  p.40)  noted  that  a 
two-dimensional  model  would  be  necessary  for  temperature 
prediction  of  bridges  at  other  than  an  east-west  orientation. 

They  developed  a  finite-element  model  for  their  analysis,  but  for 
the  purposes  of  this  study,  a  finite-element  analysis  was  thought 
to  be  too  Intensive. 

Application  of  the  TDA  models  in  the  field,  very  often  on 
hand-held  computers  (Wachtmann  et  al,  1985,  p.l5),  would  preclude 
use  of  a  model  requiring  large  amounts  of  memory  or  program 


'ii 


instructions  as  a  finite-element  model  would*  However,  it  is 


important  to  know  the  effect  of  solar  side-heating  and  determine 


if  the  one-dimensional  model  adequately  predicts  the  surface 


temperature. 


For  the  purposes  of  this  study,  it  was  decided  to  develop  a 


simple  finite-difference  two-dimensional  model  to  study  the 


effects  of  solar  side-heating  and  for  application  to  a  TDA  model 


if  appropriate.  This  development  involves  simply  applying  the 


finite-difference  technique  to  two  coordinates  of  the  Fourier 


equation,  which  gives 


Kh 


+  -:;(T 


2'‘'-i,j+l  ■  ^■^ij  '‘^l.j-l^ 


where  dx»  dy>  and  the  T^.s  are  as  shown  in  Figure  6.  If  this 


equation  is  developed  as  the  one-dimensional  model,  for  the 


temperatures  on  the  far  left  of  the  slab,  the  following  equations 


are  obtained: 


T'll  -  Tn  +  4h-2(T2i  -  Tji)  +  2hW^(Tii) 


+  4h-^(T^2  ■  ^iP  2hWy(TiP 


l-VCV.* 


T21  -  T21  +  h-^(T3i  ■  3T21  +  2Tii) 

+  -  Tjj)  +2hW  (I^j) 

dy 


t;,  -  t,,  +  h-5(T. 


•ji  “  *vi  ^  — VA  "•  T  +  T  ^ 

i+1,1  ii  i-i,r 

*  *''^<^2  -  hl^  * 

'n-1.1  -  Vl.l  *  -  3Tn.i,i  +  V2_i) 

*  ‘'’^'Vi.2  -  Vl,l^  *  2'’VVl,l> 


’'ni  ■  ’’ni  *  1  -  f»t)  +  2hW^(T^j) 


N1  Nl  ^“^'*N-1,1  "Nl'  " 

*  «»%t»2  -  ■'m)  *  ^''“ytT^i)- 

d 

y 


?»-X/  I 
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Figure  6.  Spacing  of  Temperature  Values  in  2-D  Model. 


Calculation  of  the  next  column  of  temperatures  to  the  right  would 
Involve  using  a  y^component  of  the  equation  similar  to  the 
x-component  of  the  second  equation  above*  The  y-component  would 
continue  to  change  similarly  for  each  succeeding  column  to  the 


III.  SENSITIVITY  ANALYSIS 


In  order  to  properly  apply  the  modelling  techniques  to  the 
physical  problem,  it  was  necessary  to  study  the  Impact  of  various 
factors.  Factors  such  as  the  sensitivity  of  the  models  to  changes 
in  the  finite  time  step  and  spatial  step,  key  physical  properties, 
and  certain  environmental  interactions  were  investigated.  Also, 
the  impact  of  the  two  different  approximation  techniques  was 
considered.  These  sensitivity  analyses  were  performed  using  the 
one-dimensional  roadway  structure  model. 


Physical  properties  used  in  performing  these  analyses  were 
obtained  from  the  work  by  Thepchatri  and  coworkers  (1977,  p.60). 
These  are 


Absorptivity  of  Surface  to  Solar  Radiation 
Emlssivity 


Density 


Thermal  Conductivity 


Specific  Heat 


150  Ibm/ft’ 
2403  kg/m* 


0.81  Btu/ft/hr/°F 
1.4  W/m/°C 

0.23  Btu/lbm/°F 
963  W-s/kg/°C 


Generally,  these  models  were  set  up  as  If  reflected  radiation 
from  the  water  was  incident  on  the  lower  surface  of  the  bridge# 
When  emitted  radiation  was  also  considered,  the  water  was  assumed 


to  be  at  a  constant  temperature  throughout  the  day.  Based  on 
Informal  studies  (Marlow,  1986)  this  temperature  was  chosen  to  be 
220C.  The  reflectivity  and  emissivity  of  the  water  were  assumed 
to  be  0.05  and  0.95  respectively  (Wolfe  and  Zissls,  1978,  p.3-106)< 


It  was  also  necessary  to  choose  an  initial  temperature 
profile  for  the  bridge.  Based  on  previous  studies  (Thepchatrl  et  al, 
1977,  p.60),  the  initial  temperature  profile  of  the  bridge  was 
assumed  to  be  constant  and  at  the  ambient  air  temperature  at  the 
start  time.  This  initial  temperature  was  19.25°C  at  the  starting 
time  of  the  prediction,  6:07  AM. 


Time  and  Spatial  Sensitivity: 

In  any  finite~difference  technique,  the  size  of  the  finite 
steps  is  Important.  Too  large  a  step-size  results  in  unacceptable 
inaccuracies  while  too  small  of  a  step  may  require  an  excessive 
amount  of  computer  run  time.  The  key  is  to  find  the  largest 
step-size  which  gives  an  acceptable  degree  of  accuracy. 


Another  consideration  with  the  models  developed  here  is 
choosing  a  combination  of  time  and  spatial  step-sizes  which  will 
not  result  in  numerical  instabilities.  In  this  case,  numerical 
Instabilities  will  result  when  the  quantity 


d 


m 


becomes  less  than  zero.  Therefore  to  avoid  numerical 
Instabilities  the  time  step  and  layer  thickness  should  be  selected 
such  that 


To  Investigate  the  sensitivity  of  the  models  to  variations  in 
time-step,  a  three  layer  model  of  the  roadway  structure  was  used. 
Since  the  roadway  structure  is  7.5  inches  thick  (see  schematics  in 
Appendix  B),  this  gives  a  layer  thickness  of  2.5  inches,  or  0.0635 
meters.  Time-steps  of  one,  five,  ten,  and  fifteen  minutes  were 
compared. 

The  results  of  the  time-step  sensitivity  analysis  are  given 
in  Table  1,  which  displays  the  results  in  increments  of  one  hour 
after  the  initial  time.  This  shows  very  little  change  as  the 
time-step  is  varied  from  one  to  15  minutes,  with  the  variation  in 
temperature  always  being  less  than  two  degrees  Celsius  and  most 
of  the  time  less  than  one  degree.  It  is  interesting  to  note, 
however,  that  the  maximum  difference  in  the  magnitude  of  the  peak 
temperatures  is  less  than  0.6*^C. 

Table  1  also  shows  consistent  results  between  the  Runge-Kutta 
and  the  finite-difference  models,  with  differences  at  a  given  time 
and  time-step  for  the  most  part  being  less  than  one  degree 
Celsius  (with  one  notable  exception  at  19:07  using  a  15  minute 
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Clme  step  where  the  difference  Is  1.85^C>.  For  Che  five  alnute 
and  one  alnute  doe-steps,  Che  difference  between  the  two  aodels 
Is  less  Chan  0.3°C. 

The  doe  of  occurence  of  Che  peak  teoperacure  varies  as  the 
Clae-sCep  varies  but  Is  consistent  between  the  two  aodels  except 
for  Che  15  olnute  doe  step*  The  variation  In  the  peak  occurence 
Is  soall  for  the  doe-steps  up  to  10  olnutes  and  could  be  a  result 
of  the  coarseness  of  the  aodels.  The  Cloe  of  the  peak  using  the 
15  minute  time-step  with  Che  Runge-KuCC«  model  is  also  reasonably 
close  Co  Che  other  time-steps,  but  the  15  minute, 

finite-difference  result  is  too  far  removed  from  the  ocher  results. 

To  study  the  sensitvicy  of  the  models  to  variations  in  layer 
thickness,  a  one  minute  doe-step  was  used  with  three,  five,  ten, 
and  fifteen  layer  models.  Again,  both  Che  Runge-KuCCa  and 
finite-difference  models  were  Investigated.  Layer  thicknesses 
corresponding  to  the  number  of  layers  used  are  0.0636  m,  0.0381  m, 
0.0190  m,  and  0.0127  m  respectively.  It  should  also  be  noted  Chat 
15  layers  is  the  most  that  can  be  used  with  a  one  minute  time-step 
and  avoid  numerical  .  istabillcles  with  this  model  and  the  given 
physical  properties.  A  smaller  time  step  could  be  used  with  more 
layers,  but  one  minute  was  Che  minimum  used  in  this  study  because 
the  weather  data  avallabe  was  given  at  approximately  one  minute 


intervals 


Table  2  shows  the  results  of  the  spatial  sensitivity  study . 
These  results  indicate  that  varying  the  number  of  layers  produces 
a  greater  variation  in  the  temperature  prediction  than  varying  the 
time-step.  In  this  case*  a  variation  of  as  much  as  four  degrees 
Celsius  Is  seen  occurlng  at  13:07.  However,  the  variation  1$ 
less  than  two  degrees  if  the  three  layer  model  is  not  considered. 
Also,  the  variation  In  the  magnitude  of  the  peak  temperature  is  as 
much  as  2.75^C,  but  only  l.4°C  if  the  three  layer  model  is  not 
considered.  If  just  the  ten  and  15  layer  models  are  considered, 
the  variation  Is  less  than  O.S^C.  It  is  also  noted  that  the  time 
of  occurence  of  the  peak  is  consistent  with  all  the  models  used  in 
the  spatial  sensitivity  study. 

The  consistency  between  the  finite-difference  and  the 
Runge-Kutta  models  is  better  as  the  number  of  layers  is  varied 
than  if  the  time-step  is  varied.  The  data  from  Table  2  indicates 
that  for  a  given  number  of  layers  the  variation  between  the  two 
models  at  any  time  is  less  than  0.12°C. 

The  results  of  the  time  and  spatial  sensitivity  studies 
indicate  that  reasonable  accuracy  can  be  obtained  with  either  the 
finite-difference  or  the  Runge-Kutta  model.  However,  to  maintain 
the  accuracy,  the  time-step  should  be  less  than  five  minutes  and 
the  number  of  layers  should  be  greater  than  five.  For  further 
sensitivity  studies,  only  the  finite-difference  model  was  used. 


Senslclvlcy  Co  Changing  Dlffuslvlc 


Concrete  bridges  such  as  the  one  considered  here  typically 
have  steel  reinforcements.  The  amount  of  steel  used  In  any 
particular  bridge  may  vary,  which  would  greatly  alter  the 
dlffuslvlty  factor  used  In  Che  model.  This  results  because  the 
thermal  conductivity  of  steel  Is  on  the  order  of  30  times  greater 
chan  chat  of  concrete,  while  the  product  of  the  density  and 
specific  heat  of  steel  Is  only  about  1.6  times  greater  chan  chat 
of  concrete. 


For  example,  average  values  of  concrete  thermal  properties 
applicab'i.e  to  steel  reinforced  bridges  are  (Thepchatrl  et  al. 


1977,  p.60) 


Thermal  conductivity 


Density 


Specific  heat 


0.81  Btu/ft/hr/  F 

1.4  W/mrc 

150  Ib/ft^ 

2403  kg/m^ 

0.23  Btu/lb/°F 
963  W-s/kg/  C 


Comparable  values  for  steel  are  (Rohsenow  and  Hartnett,  1973, 
pp. 2-67, 68) 


Thermal  conductivity 


Density 


Specific  heat 


26.6  Btu/ft/hr/°F 
46  W/m/°C 

500  Ib/ft^ 

8009  kg/m^ 


0.11  Btu/lb/°F 
460  W-s/kg/°C 


As  used  In  the  model »  the  surface  properties  of  the  bridge 
are  all  divided  by  the  product  of  the  density  and  specific  heat* 


a'B  a*#  MP*  M»* ~^jm~ jnj  m^^rn.  ar^ _ j~i 


.  BVa  g 


However,  this  product  is  very  large  for  both  coacrete  and  steel 
while  the  surface  properties  are  relatively  very  soiail  numbers. 
From  this  It  was  concluded  that  the  effect  of  varying  amounts  of 
steel  In  the  bridge  struture  would  be  almost  entirely  from  the 
dlffuslvlty  term. 


Therefore,  to  study  this  effect,  the  model  was  run  for  cases 
where  the  dlffuslvlty  term  Is  multiplied  by  factors  of  two,  five, 
and  eight.  A  five  layer,  one  minute  time-step  model  was  used. 

The  one  minute  time-step  was  selected  to  avoid  numerical 
Instabilities  with  the  larger  dlffuslvlty  terms. 


Table  3  shows  the  results  of  this  analysis.  The  magnitude 
and  time  of  occurence  of  the  peak  temperature  Is  of  most 
significance  here.  It  can  be  seen  that  the  peak  temperature  Is 

reduced  by  almost  nine  degrees  Celsius  when  the  dlffuslvlty  is 
eight  times  Its  normal  value,  and  the  time  of  the  peak  Is  delayed 
until  almost  1.3  hours  later.  Also,  the  results  using  a  factor  of 
five  are  very  close  to  those  obtained  using  a  factor  of  eight, 
Indicating  a  levellng-out  of  this  effect.  A  slight  lowering  of 
the  peak  temperature  without  shifting  the  time  of  occurence  Is 
obtained  using  a  factor  of  two  times  the  dlffuslvlty. 


The  effect  of  changing  dlffuslvlty  Is  shown  graphically  in 
Figure  7  where  the  resulting  temperatures  for  cases  of  one  times 


TEHPERATURE  (°C) 


TABLE  3 

SENSITIVITY  TO  CHANGING  DIFFUSIVITY 


DIFFUSIVITY  -  TEMPERATURES  -  “c 


FACTOR 

IX 

2X 

5X 

8X 

TIME 

7:07 

17.64 

17.90 

18.26 

18.41 

8:07 

23.40 

22.57 

21.69 

21.30 

9:07 

30.29 

28.30 

26.12 

25.26 

10:07 

35.63 

32.93 

30.02 

28.98 

11:07 

43.31 

39.44 

35.42 

34.01 

12:07 

51.06 

46.13 

41.14 

39.44 

13:07 

54.57 

49.72 

44.90 

43.34 

U:07 

56.84 

52.30 

47.89 

46.48 

15:07 

56.96 

53.10 

49.37 

48.21 

16:07 

55.79 

52.66 

49.76 

48.88 

17:07 

51.26 

49.49 

47.92 

47.51 

18:07 

44.05 

43.85 

43.96 

44.10 

19:07 

34.23 

35.86 

37.68 

38.40 

20:07 

26.67 

29.33 

32.17 

33.21 

21:07 

22. '22 

25.14 

28.08 

29.06 

22:07 

18.85 

21.76 

24.57 

25.46 

23:07 

16.29 

19.01 

21.56 

22.37 

MIDNIGHT 

14.63 

17.10 

19.39 

20.14 

TIME  OF 

PEAR  TEMP. 

14:41 

14:41 

16:02 

16:08 

PEAK  TEMP. 

57.68 

53.39 

49.80 

48.91 

NOTE;  TEMPERATURES  MODELLED  USING  I  MINUTE  TIME- 
STEP,  5-LAYER,  FINITE-DIFFERENCE  MODEL. 


Effect  of  Changing  Diffusivity  on  Surface  Temperatures 
with  Reflected  Radiation  only  Incident  on  Lower  Surface 


and  eight  times  the  dlffuslvlty  are  plotted  against  time.  Both 
the  top  surface  and  lower  surface  temperatures  are  plotted  to  show 
how  they  approach  each  other  more  closely  with  a  higher 
dlffuslvlty. 

The  shifting  of  the  peak  Is  likely  due  to  the  effect  of  the 
radiation  Incident  on  the  lower  surface.  This  Incident  radiation 
serves  to  heat  the  lower  surface,  and  with  the  higher  dlffuslvlty, 
this  heat  Is  more  readily  conducted  to  the  top  surface.  The  model 
used  for  this  sensitivity  analysis  considered  only  reflected 
radiation  from  the  water  onto  the  lower  surface.  Varying  effects 
of  radiation  Incident  on  the  lower  surface  were  also  considered. 

Lower  Boundary  Sensitivity; 

Variations  In  the  surroundings  or  structure  of  a  particular 
bridge  may  affect  the  amount  of  radiation  received  by  Its  lower 
surface.  In  the  sensitivity  analyses  discussed  so  far,  the  models 
considered  only  reflected  radiation  from  the  water  Incident  on  the 
lower  surface.  However,  radiation  emitted  by  the  water  may 
greatly  Impact  the  radiation  heat  transfer  to  the  lower  surface, 
or  the  structure  may  be  such  that  no  radiation  from  the 
surroundings  Is  Incident  on  the  lower  surface. 


For  most  TDA  applications,  the  temperature  of  the  lower 
surface  Is  of  no  consequence  except,  as  shown  In  the  previous 


subsectloa,  for  how  It  affects  the  top  surface  temperature* 
Therefore,  it  is  important  to  know  how  varying  amounts  of 
radiation  from  the  surroundings  incident  on  the  lower  surface  will 
affect  the  top  surface  temperature. 

To  study  this  effect,  the  finite-difference  model  was  run  for 
cases  of  no  radiation,  reflected  radiation  only,  emitted  radiation 
only,  and  both  reflected  and  emitted  radiation  incident  on  the 
lower  surface*  To  see  the  relation  of  the  top  surface  temperature 
with  that  of  the  lower  surface,  the  results  of  these  cases  are 
shown  in  graphical  form  in  Figures  7-10  where  both  the  top  and 
lower  surface  temperatures  are  plotted  for  each  case* 


LOCAL  TIME  (hrs.) 


Figure  8*  Effect  of  No  Radioslty  Incident  on  Lower  Surface. 


It  Is  noted  from  these  figures  that  large  Increases  are  seen 
in  the  lower  surface  temperatures  (as  much  as  20^C  for  the  extreme 
cases)*  However,  the  Increase  in  the  top  surface  temperature  Is 
only  about  four  degrees  Celsius. 
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IV.  MODELLING  RESULTS 


An  appropriate  form  of  the  one-dlmensional  finite-difference 
model  was  selected  to  make  a  prediction  of  the  bridge  surface 
temperatures  for  the  day  of  the  measurements.  The  selection  of 
this  model  was  based  on  the  results  of  the  sensitivity  analyses 
and  estimates  of  the  surface  temperatures  made  from  the  radiance 
measurements.  Using  this  prediction  as  a  baseline,  a 
two-dimensional  finite-difference  model  was  also  used  to  predict 
the  surface  temperature  distribution  for  comparison  to  the 
one-dimensional  results.  The  two-dimensional  model  was  also  used 
to  Investigate  the  side-heating  impact  of  a  north-south  oriented 
bridge.  Sample  program  listings  for  both  the  one-dimensional  and 
two-dimensional  models  can  be  found  in  Appendix  C. 

Just  one  of  the  developed  one-dimensional  models  was  carried 
through  this  process  based  on  the  results  of  the  sensitivity 
analyses.  These  showed  that  the  variation  between  results  of  the 
two  models  was  insignificant  provided  the  proper  time  and  spatial 
spacing  was  used.  The  finite-difference  model  was  selected 
because  it  is  more  straight-forward  and  easier  to  Implement.  This 
model  was  used  with  a  five  minute  time-step  and  five  layers,  and 
again  the  starting  time  was  6:07  AM  and  the  Initial  temperature 
distribution  was  assumed  to  be  a  constant  19.25^C. 


Estimates  of  Actual  Temperatures; 

The  radiance  measurements  made  by  AFATL  on  17  May  1984 
(TABILS  I)  provided  a  means  to  estimate  the  actual  temperatures  of 
the  bridge.  Scenes  of  the  bridge  in  the  7,0  to  14,0  micrometer 
band  taken  at  various  times  during  the  day  were  analyzed  using  the 
process  described  by  Goldman  and  coworkers  (1983,  pp, 50-53).  This 
process  provided  a  value  for  the  radiance  at  the  detector,  Iq. 
Using  the  weather  data  taken  on  the  day  of  measurements  and  the 
system  spectral  response,  a  transmittance  value,  t,  combining  the 
path  transmittance  and  system  spectral  response  was  found.  Also, 
a  value  for  the  radiance  of  the  path.  Ip,  was  found  using  the 
weather  data. 

With  an  estimate  of  these  values,  the  radiance  at  the  source, 
Ig,  could  be  estimated.  This  is  given  by 


t 


However,  this  estimate  of  the  source  radiance  Includes  both 
emitted  and  reflected  radiation.  To  estimate  the  source 
temperature,  the  value  of  the  emitted  radiation  alone  needs  to  be 
known. 


Goldman  and  coworkers  (1983,  pp. 68-77)  note  that  the 
contribution  of  reflected  radiation  can  be  significant.  However, 
their  data  Indicates  that  for  an  object  with  an  eraisslvity  of  0,9 


in  the  temperature  range  considered  here,  reflected  radiation 
accounts  for  approximately  5Z  of  the  source  radiation.  So,  the 
value  found  for  the  source  radiation  using  the  above  equation 
reduced  by  5%  was  used  as  the  value  for  the  emitted  radiation,  Igs* 

The  temperature  estimate  based  on  the  value  l£3  was  made 
using  an  iterative  process.  First,  a  temperature  was  assumed  from 
which  the  percentage  of  radiation  in  the  given  band  was  found. 
Using  this  value,  an  emissivlty  value  of  0.9,  the  value  of  Igg, 
and  the  Stefan-Boltzmann  Law;  a  new  temperature  was  found.  This 
process  was  repeated  until  the  temperatures  agreed.  The  results 
of  this  analysis  are  given  in  Table  4. 

It  should  be  noted  here  that  these  temperature  esclmates  can 
not  be  regarded  as  true,  measured  values.  No  reliable  temperature 
measurements  of  the  bridge  surface  on  the  day  of  the  radiance 
measurements  were  available  for  this  study.  For  this  reason,  the 
radiance  values  were  used  to  identify  trends  in  the  daily  cycle  of 
the  bridge  surface  temperature.  Identifying  these  trends  was 
useful  in  determining  appropriate  environmental  interactions  to 
consider. 

One  Dimensional  Model  Results; 

The  finite-difference  model  using  a  five  minute  time-step  and 


five  layers  was  run  for  both  the  roadway  and  the  piling.  The 


TABLE  4 

ESTIMATED  EMITTED  RADAIANCE  AND  TEMPERATURE  OF 
THE  DESTIN  BRIDGE  FROM  17  MAY  1984  DATA 


(HH:MM:SS) 


PATH  RAD.*  TRANS. 


ROADWAY  RAD.*  &  TEMP.  PILING  RAD.*  &  TEMP. 
APP.  SOURCE#  EST  TEMP  APP.  SOURCE#  EST  TEMP 


thicknesses  of  the  roadway  and  the  piling  are  7.5  Inches  and  24 
Inches  respectively  (see  Appendix  B).  These  models  were  run 
considering  no  reflected  or  emitted  radiation  from  the  water 
Incident  on  the  bridge  structure.  The  reason  for  this  approach  Is 
because  of  the  trends  Identified  from  the  estimated  temperatures 
and  the  lower  boundary  sensitivity  analysis. 

The  results  of  the  model  chosen  are  plotted  In  Figure  11, 
which  also  shows  the  estimated  temperature  data.  This  shows  that 
the  model  predicts  the  piling  becomes  warmer  than  the  roadway  top 
surface  at  about  the  same  time  as  Indicated  by  the  estimated 


LOCAL  TIME  (hrs.) 


Figure  11.  Roadway  and  Piling  Surface  Temperature  Predictions 

(No  Radloslty  Incident  on  Lower  Surface)  Compared  to 
Temperature  Estimates. 


temperature  data.  This  was  true  even  when  emitted  and  reflected 


Figure  13,  Roadway  and  Piling  Surface  Temperature  Prediction 
with  Emitted  and  Reflected  Radiation  Incident  on 
Lower  Surface. 


As  noted  in  the  sensitivity  analysis,  there  is  little 
difference  in  the  results  for  the  top  surface  when  both  emitted 
and  reflected  radiation  are  assumed  incident  on  the  lower 
surface.  Adding  the  effect  of  emitted  and  reflected  radiation 
drives  the  temperature  of  the  piling  nearly  15°C  higher  in  the 


late  afternoon. 


Two-Dimensional  Model  Results: 


A  two-dimensional,  finite-difference  model  was  developed  In 
order  to  study  the  effects  of  side  heating  In  both  the  east-west 
orientation  of  the  Destln  bridge  and  In  a  supposed  north-south 
orientation.  Both  the  roadway  and  the  piling  were  modelled  In 
two-dlmenslons ,  and  a  five  minute  time-step,  as  used  with  the 
one-dlmenslonal  model,  was  selected.  The  roadway  Is  48  ft.  and 
7.5  In.  wide  and  the  pilings  are  2  ft.  square  (see  Appendix  B). 
Spacing  for  the  roadway  was  five  layers  across  the  thickness  like 
the  one-dlmenslonal  model  and  15  layers  across  the  width.  For  the 
piling,  five  layers  was  again  used  across  the  thickness  and  the 
width  was  divided  Into  10  layers.  This  spacing  produced  cells  of 
0.0381  by  0.988  meters  In  the  roadway  and  cells  of  0.122  by  0.061 
meters  In  the  piling. 

When  modelled  In  an  east-west  direction,  the  two-dimensional 
roadway  model  shows  that  there  Is  very  little  variation  In  the 
surface  temperature  across  the  width  of  the  roadway.  Table  5 
shows  the  results  of  this  model  compared  to  the  results  of  the 
one-dlmenslonal  model.  From  this  It  Is  seen  that  the  temperature 
of  the  Interior  layers  across  the  surface  matches  that  predicted 
by  the  one-dlmenslonal  model.  Only  In  the  end  layers  and  end 
surfaces  Is  a  slight  difference  seen. 
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As  shown  in  Table  6,  more  variation  is  seen  in  the  surface 
temperature  of  the  piling  when  it  is  modelled  in  two  dimensions* 
This  was,  of  course,  expected  since  the  piling  is  a  much  smaller 


structure  than  the  roadway.  The  center  temperature,  however, 
stays  very  close  to  that  predicted  by  the  one-dimensional  model. 
Variations  in  temperature  from  the  center  to  the  ends  are  less 
than  2.5°C.  Therefore,  an  average  temperature  across  the  surface 
of  the  piling  would  not  be  very  different  than  that  predicted  by 
the  one-dimensional  model,  indicating  that  a  one-dimensional 
prediction  could  be  used  to  evaluate  the  response  of  an  infrared 
detector  with  reasonable  accuracy. 


Since  the  Destin  Bridge  runs  east-west,  the  sides  and  the 
pilings  underneath  are  subjected  to  a  minimum  of  solar  radiation 
However,  for  a  bridge  not  oriented  in  an  east-west  direction, 
solar  side  heating  could  have  a  significant  effect  on  the  surface 
temperature.  To  investigate  the  effect  of  solar  radiation  on  the 
sides  and  the  piling,  the  two-dimensional  models  were  run  for  the 
bridge  as  if  it  were  at  a  north-south  orientation. 


The  measurements  of  the  solar  radiation  given  in  the  weather 
data  are  for  the  radiation  incident  on  a  horltzontal  surface.  For 
the  purpose  of  the  north-south  model,  it  was  necessary  to  estimate 
the  solar  radiation  incident  on  a  vertical  surface. 


TIME 

7:07 

1-D 

19.31 

00 

o 

19.53 

9:07 

20.07 

10:07 

20.58 

11:07 

21.10 

12:07 

21.73 

13:07 

22.40 

14:07 

23.00 

15:07 

23.60 

16:07 

24.05 

17:07 

24.40 

18:07 

24.25 

19:07 

23.88 

20:07  . 

23.94 

21:07 

23.81 

22:07 

23.49 

23:07 

23.09 

23:57 

22.78 

TEMPERATURES  -  C 
19.31  19.31  19.31 


19.38  19.33  19.32  19.31  19.31  19.31  19.31  19.31  19.31  19.32  19.33  19.38 
19.79  19.62  19.55  19.53  19.53  19.53  19.53  19.53  19.53  19.55  19.62  19.79 
20.77  20.36  20.14  20.09  20.07  20.07  20.07  20.07  20.09  20.14  20.36  20.77 

21.56  21.11  20.76  20.63  20.59  20.58  20.58  20.59  20.63  20.76  21.11  21.56 

22.39  21.89  21.41  21.21  21.13  21.11  21.11  21.13  21.21  21.41  21.89  22.39 

23.50  22.81  22.19  21.92  21.80  21.76  21.76  21.80  21.92  22.19  22.81  23.50 

24.50  23.76  23.03  22.67  22.51  22.45  22.45  22.51  22.67  23.03  23.76  24.50 
25.32  24.61  23.80  23.38  23.17  23.09  23.09  23.17  23.38  23.80  24.61  25.32 
26.10  25.42  24.56  24.08  23.83  23.73  23.73  23.83  24.08  24.56  25.42  26.10 
26.65  26.01  25.15  24.63  24.36  24.24  24.24  24.36  24.63  25.15  26.01  26.65 
27.02  26.45  25.62  25.09  24.78  24.65  24.65  24.78  25.09  25.62  26.45  27.02 
26.28  26.14  25.52  25.02  24.71  24.57  24.57  24.71  25.02  25.52  26.14  26.28 
25.48  25.47  25.09  24.69  24.40  24.26  24.26  24.40  24.69  25.09  25.47  25.48 

25.57  25.41  25.08  24.76  24.51  24.37  24.37  24.51  24.76  25.08  25.41  25.57 
25.19  25.14  24.90  24.62  24.41  24.29  24.29  24.41  24.62  24.90  25.14  25.19 
24.45  24.61  24.49  24.28  24.10  24.00  24.00  24.10  24.28  24.49  24.61  24.45 
23.73  23.95  23.97  23.85  23.71  23.63  23.63  23.71  23.85  23.97  23.95  23.73 
23.21  23.45  23.56  23.50  23.40  23.33  23.33  23.40  23.50  23.56  23.45  23.21 


FIRST  AND  LAST  2-D 
APPROXIMATELY  0.03 
APART. 


TEMPERATURES  ARE  AT  ENDS.  2ND 
METERS  FROM  ENDS.  ALL  OTHERS 


AND  IITH  TEMPERATURES  ARE 
ARE  APPROXIMATELY  0.06  METERS 


To  estimate  the  solar  radiation  incident  on  a  vertical 


surface,  it  was  noted  that  the  measurements  of  the  solar  radiation 
were  approximately  sinusoidal  with  respect  to  the  angle  between 
the  incoming  radiation  and  the  normal  to  the  surface,  peaking 
around  mid-day  when  the  sun  would  be  directly  overhead.  It  was 
assumed  then  that  for  a  vertical  surface  the  magnitude  of  the 
incident  solar  radiation  would  follow  a  cosine  curve  with  respect 
to  the  same  angle. 

The  sun  was  assumed  to  be  directly  overhead  at  a  time  halfway 
between  the  times  of  the  first  and  last  non-zero  solar  radiation 
measurements.  Then  the  angle  between  the  horltzontal  surface 
normal  and  the  solar  radiation,  @,  was  assumed  to  change  at  a 
constant  rate.  The  magnitude  of  the  solar  radiation  incident  on  a 
vertical  surface  was  then  estimated  using  the  following  equation: 


W  _  *  W  cos0/sin0  . 
vert  sun '  ' 


It  was  assumed  that  the  bridge  structure  would  obscure  solar 
radiation  from  the  side  of  the  piling  facing  Inward.  Therefore, 
two  north-south  piling  models  were  used:  one  with  the  piling 
facing  the  morning  sun  and  one  facing  the  afternoon  sun. 

Table  7  shows  a  comparison  of  the  one-dimensional  model 
results  to  the  results  of  modelling  the  roadway  in  a  north-south 
direction.  The  north-south  model  was  oriented  so  the  left  side 
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(corresponding  to  the  2-D  temperatures  on  the  left  of  Table  7) 
faced  the  morning  sun.  Here,  it  is  again  seen  that  the  interior 
layer  temperatures  in  the  two-dimensional  model  match  the 
temperatures  predicted  by  the  one-dimensional  model.  However,  the 
effect  of  the  solar  side  heating  causes  more  difference  between 
the  end  and  interior  temperatures.  This  difference,  though,  is 
still  not  significant  in  any  of  the  interior  layers. 

Tables  8  and  9  show  the  results  of  modelling  the  piling  as  if 
it  faced  the  morning  and  afternoon  sun  respectively.  The  center 
temperatures  are  in  close  agreement  with  those  modelled  with  a 
north-south  oriented  one-dimensional  model,  comparable  to  the 
difference  seen  for  the  east-west  orientation. 

However,  Tables  8  and  9  also  show  more  variation  between  the 
center  and  end  temperatures.  The  variation  is  not  significantly 
higher,  though,  than  that  shown  in  Table  6  for  the  east-west 
model.  Table  8  indicates  that  in  the  morning  for  a  piling  facing 
the  morning  sun  the  variation  between  the  temperatures  at  the 
center  and  at  the  ends  is  less  than  about  five  degrees  Celsius. 
This  maximum  occurs  early  at  8:07,  and  the  center  temperature  is 
higher  than  the  end  temperature  where  for  the  east-west  model  the 
end  was  warmer  than  the  center.  This  results  from  the  piling 
being  warmed  by  the  sun  in  the  morning  and  that  heat  being 
transferred  to  the  surroundings  more  readily  at  the  ends.  In  the 
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TABLE  8 

COMPARISON  OF  PILING  1-D  AND  N-S  MODEL  RESULTS 
FACING  MORNING  SUN 
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TABLE  9 

COMPARISON  OF  PILING  1-D  AND  N-S  MODEL  RESULTS 
FACING  AFTERNOON  SUN 


-  TEMPERATURES  -  °C 
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east-west  model  the  surroundings  are  warmer  than  the  piling  so 
heat  is  transferred  more  readily  to  the  ends.  This  difference 
decreases  slowly  until  the  ends  become  warmer  than  the  piling  at 
around  16:07. 


Table  9  shows  a  slightly  different  scenario.  The  results  for 
the  north-south,  afternoon  sun  facing  model  naturally  follows  the 
results  of  the  east-west  model  in  the  morning.  At  13:07,  when  the 
sun  is  starting  to  shine  on  the  piling,  the  results  start  to 
become  different.  The  center  temperature  starts  to  approach  that 
of  the  end  and  is  warmer  by  17:07.  After  that,  the  cooling  of  the 
night-time  is  much  like  that  for  the  east-west  model. 


One  final  comparison  that  can  be  made  is  to  see  if  in  a 
north-south  orientation  the  piling  would  become  warmer  than  the 
roadway,  thereby  being  more  of  a  target.  Comparing  Tables  7  and 
8,  it  is  seen  that  a  piling  facing  the  morning  sun  would  be  much 
warmer  than  the  roadway  until  about  mid-morning.  Table  9  shows 
that  a  piling  facing  the  afternoon  sun  would  not  become  warmer 
than  the  roadway  until  very  late  in  the  afternoon,  but  it  will 
stay  much  warmer  than  the  roadway  until  late  in  the  night.  This 
la  partly  a  result  of  the  roadway  losing  much  of  its  heat  to  the 
night  sky  while  the  piling  is  not  subjected  to  this  phenomena 
since  it  does  not  directly  view  the  night  sky. 


It  should  be  noted  from  Tables  5  and  6  that  even  when  the 
piling  is  subjected  to  no  solar  radiation  during  the  day  it  is 
warmer  than  the  roadway  late  at  night*  However,  being  subjected 
to  an  afternoon  sun  causes  the  piling  to  become  warmer  than  the 
roadway  much  sooner* 


V,  CONCLUSIONS  AND  REMARKS 


This  study  has  examined  the  numerical  application  of  the 
Fourier  heat  equation  to  the  problem  of  surface  temperature 
prediction  of  a  bridge-structure.  The  reason  for  approaching  this 
problem  was  to  provide  a  means  for  estimating  the  apparent 
contrast  of  a  bridge  using  a  Tactical  Decision  Aide  model. 

From  the  results  of  this  study,  several  conclusions  can  be 
drawn.  These  are  as  follows: 

1.  Either  a  Runge-Kutta  or  a  full  finite-difference 
approximation  provides  nearly  the  same  prediction  of  surface 
temperature. 

2.  A  one-dimensional  model  provides  nearly  the  same 
prediction  as  would  a  two-dimensional  model. 

3.  The  top  surface  temperature  of  a  roadway  structure  is 
affected  very  little  by  the  conditions  of  the  lower  surface. 

4.  The  pilings  of  a  bridge  oriented  in  a  north-south 
direction  will  provide  a  more  visible  target  than  the  roadway  in 
either  the  early  morning  or  late  afternoon,  depending  on  what 
direction  the  piling  faces. 

5.  Prediction  of  the  roadway  surface  temperature  is 
negligibly  affected  by  orientation. 
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It  was  not  possible  In  this  study  to  determine  the  degree  of 
accuracy  to  which  the  models  predicted  the  surface  temperatures. 
Trends  in  the  variation  of  the  surface  temperature  throughtout  the 
day  were  approximately  matched  based  on  radiance  measurements  of  a 
particular  bridge.  No  reliable  ground-truth  data  was  available  to 
compare  to  the  predictions.  More  extensive  measurements  of  a 
particular  bridge  or  bridges,  including  both  radiance  measurements 
and  temperature  measurements,  would  be  necessary  to  provide  enough 
data  to  determine  the  accuracy  of  the  predictions. 

This  study  only  considered  the  interaction  of  the  bridge  with 
a  particular  background — water.  The  effects  of  other  backgrounds 
such  as  vegatation,  rock,  bare  earth  or  various  combinations 
should  also  be  considered. 

The  use  of  different  materials  may  also  have  an  impact  on  the 
surface  temperature  prediction.  The  bridge  that  was  the  focus  of 
this  study  was  mostly  concrete,  and  some  consideration  was  given 
to  the  effect  of  steel  reinforcements.  Consideration  should  be 
given  to  a  bridge  which  is  mostly  or  all  steel.  Also,  a  wooden 
bridge  should  be  considered  from  the  standpoint  that  it  may  be  a 
military  target. 
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A.  THERMOGRAPH  SCENES  OF  THE  DESTIN  BRIDGE 


Figure  A-1  shown  here  contains  thermograph  scenes  of  the 
Destln  Bridge  taken  on  17  May  1984  with  an  AGA  680  Thermovision 
Camera*  These  scenes  are  In  the  7.0  to  14*0  micrometer  band* 


a*  Scene  at  5; 27: 40* 


b*  Scene  at  5:46:24. 


c*  Scene  at  6:46:36* 


d*  Scene  at  9:13:0L 


Figure  A-1*  Thermographs  Scenes  of  the  Destln  Bridge* 
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B.  EXCERPTS  FROM  DRAWINGS  OF  DESTIN  BRIDGE 


Figures  B-1  and  B-2  are  sections  from  engineering  drawings  of 
the  Destin  Bridge  showing  dimensions  of  parts  of  interest. 


Figure  B-1.  Drawing  Section  Showing  Piling  Dimensions 
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C.  SAMPLE  PROGRAM  LISTINGS 


»  1  Vin  TIME  STEP 
»  5-LAYEF  MODEL 

*  FIMITF  CIFFEFENCE  APFPC XIMAT ION 
»  REFLECTED  &  EMITTED  RADIATION  ON  LCViEF  SURFACE 

DIMENSION  TS(I074) rT(5,1074),TL(X074) 

H=60. 

N  =  5 

JENDS1074 

ALFHA1s2.16410.*4C-7.0) 

ALPHA2s3,89»10.»»(-7,C) 

ALFHA3sl.63«10.<4C-6.0) 

ALFHA45:7.30410,^<(-7,0  3 
ALFHA5s6.05410.*<(-7,05 
ALFHA6al,95*10.*< (-S.O) 

EW=l.S9<10.4<(-4.0) 

SG=5.6697*10,*<(-8,0) 

TMs2i998 

IcilNT(Cl.+FLOAT(N))/2.) 

C=O.62S/3.201 

D=C/FLOAT(N) 

AlsALPHAl/D 

A2=ALPHA2/C 

A3SALPHA3/C 

A4=ALPHA4/C 

A5sALPHA5/D4*2. 

A6SALPHA6/C 

EVisEU/0 

TS(1)  =  19. 25490 
DC  5  1=1. N 

5  TCI, n  =  l9. 25490 
TL(l)sl9.25490 

OPEN  (UNIT=9. FILE= 'EWE AT HR, DAT  *, FORMS •FORMATTED  *. ST ATUS=* CL C  * ) 
OPEN (UNITsl 6, FILES •TEMFPAD2. DAT FORMS 'FORMATTED  * ,STATUS='NEW ' ) 
TMINsH/60. 

WPITE(10,597)TMIN.N 

597  FCRMATCX, 'FDIFF  APFRX , ' , F5 . 1 , ’  MIN  TS,',I3,'  LAYERS,  FULL  RAC'/) 
WRITE(10,598) 

598  FCRMATCX, 'TIME' ,3X, 'TEMPERATURES  IN  DEGREES  C) 

WP2TEC10,599) 

599  FORMAT C2X, ' SEC ' , SX , ' TS ' , 1 0 X , < TC ' , 1 OX , ' XL ' ) 

DC  10  Isl.JEND-I 

WRITE(10,60C)TH.T£(I) ,TCIC,I) ,TL(I) 

600  FCRMATCX. F6.0,5F12,5) 

99  READC9.400)TI:vA.TA,WN,WY 
400  FCRMAT(F7.6,Fl2,3.3F15.6) 

IFCTl.EO.TM)  GO  TC  1 
GC  TC  99 

1  WSsAl»WM-FA2<CWY-SG»CTSCT)4273,15  )»»4,)  +  CA3^A4*VA)»CTA-TSCn  ) 
WLsCA3  +  A4»VA)»CTA-TLCl))  +  A2^CSG4C(TA*^273.15)»*4. 
1-CTLCI)4273.15)»*4,))4A6»C'^F  +  WY)+FW 
TSCI+l)=TS(I)+4.»A5*H»CT(l .I)-TSCI))+2.«H*WS 
TClfcl  +  nsTCI  ,I)+A54H*CTC2,I)-3.*TC1,I)+2.»TSCI)) 

CO  6  J=2,N-1 

6  TCJ,I-fl)=TCJ,I)fA5*H^(T(j4l  ,I)-2.<TCJ,I)'^T(J-1 ,1)) 
ICN,I41)sTCN,I)+A5*H»C2,»TL(I)-3.»T(N,1)4T(N-1,I)) 
TL(l4l)sTLCl542,*HtW  44.*A5tH*(TCN,I)-TLCI)) 

10  TM5TM4H 

WRITE ( 1 0. 600 ) TH,TS CJINC) ,TCIC,JEND) ,TLCJENC) 

CLCSECUN1Ts9) 

CLCSECUNITslO) 

STOP 

END 


Figure  C-1.  Program  Listing  for  1-D  Finite-Difference  Model, 


p»TSCI)) 


VO  OOO  vVr 


>•<  V  v  v> 


I*  .ta* Jb*>>  i*  «<**>»**  A*aH  Vt*  >  *  i*  A  *  A*  A  V  l«*.lb*^*J 


»  15  yv:  TIVC  CAZP 

♦  3-LAYEr  V'CDEl 

♦  '^U‘:$E-:<rTTA  Arrpf  vTr.*;*!  ir*: 

HEAL  K,K5,KL 

nTVF::;.*^lC:;  K.'!  (  n  ,»<  (  3  ,  A)  ,  KL  cn  ,YS  (7^  )  ,  Y  C3 , 7?  )  ,  YL  (72  )  ,1  O  J 

H=9CC. 

r.  =  3 

JE:ir=72 

ALFHA1=2.  J6»aO.*M-7,^) 

ALFHA2=3.39tlO,*<(-"'.C) 

AT  rVA3  =  l  C-f  .C  ) 

ALFWA4  =  7,3C'tlO.*’*  (-•'.r  1 
AL?HA5  =  ^).0t:4in,»4  (-•’.r ) 

ALFHA6  =  1.'55<1C.»*  (-?.r) 

SCs5,66'57»K. ♦♦(-?. r) 

Ty=21V9^, 

IC=ir;T(Cl,4ri.rA^('’n/2.) 

C=C. 625/3. 281 
CsC/FLOATCO 
AlsALTMAl/n 
A2  =  AL?:n2/r: 

A3=ALFHA3/n 
A4  =  AL?i:Al/r 
A5=ALFHA5/r»^2. 

AcrALTHA^^/r 
TS=19.2549C 
V5(n=T- 
nc  5  1=1. V 
1  (Dsn.ls-iso 
5  Y(I,1)=T(I) 

TL=19.2349r 
YL(n  =  TL 

CFE’)(nTlIT  =  S  ,FIL£=  *E”Er7PP.E;.T'  ^FCRI-^s  *  FCrVAlTtr  *  ,E7A1  US=  *CIC  ’  ) 
CFEN(UNIT=14  ,FILE  =  *FTF^'F2.CAT  *  ,FCPM=  'FCHFAl’ltC  *  ,£lAaUS=  ’NtW  '  ) 
7vi:^.H/6C, 

WFITE(n,597)TMIN 

597  FCF^^ATCX  , 'R-K  APF^X  ,  *  ,  F  5 . 1  ,  *  fT*'  TS,*,I3,*  LAYERS'/) 
WR1TEC10,598) 

598  FTR^^ATCX,  *TTFE*  ,SX, 'TF^^FER^TURES  IN’  CrGFEEF  C') 

WEITEC10,5q9) 

599  FCRr-*ATC7X,  'SEC*  ,Sr  ,  '  TS  ’  ,  1  '  TC  '  ,  1  OX  ,  '  Tl  *  ) 

CC  20  I=1.0END-1 

WFITEC1'>,6C':)TF,YS(I),YCIC/ID  ,yL  CT) 

600  FCR‘'AT(X.F€.0,3F12.'^) 

99  REAr(9,4C0)TI  .VA.Tm.N::  ,V:Y 
400  FCFM.ATCF7.C,Fi2.3,3F15.6) 

IFCTI.EC.TrO  GC  ir  I 
GC  TC  99 
1  X  =  2. 

CC  10  L=1.4 

VS  =  Al*W'l  +  A2^(ViY-SC4CTC4773.15)»»4.}+(A34A44VA)4ClA-l£) 

WL=(A3  +  A4»VA)*(TA-aL)'>A2»C^^G*(CTA«^27  3.15)**4. 

1-(TL4  271,15)»'»4.))^.’\6^  (V|:  +  ^Y) 

K5(L)=CU*A5*(T(1)^TS)42.*VS)*H 
K(l,l)  =  CA54(TC2)-3.’»T(n  +  2.<TS))»S 
CC  15  J=2,N-1 

15  k(j,l}  =  (:a5*ct(J+i ^7 (0 

K  CN,L)  =  (A5*(2.»TL-3.*7  (N-1)  )  )»H 

KL(t  )  =  (2.»v;L*  I.4A54  fTtrO-TT  ))»H 
1F(I  ,LT.3)  CC  TP  3 
X  =  i. 

3  TSsYSCD+KSCD/X 
CC  IG  V=i.r’ 

16  T(V)=Y(’',t)4K(M.L  )/r 
1C  TI=YL(I)+FLCL)/X 

YSCl4l)  =  yS(I)4(KS(l}4:.=>KSt2)+2.<KS(3)4K5(4))/6. 

CC  17  J=1 

17  Y(J,l4l)  =  Y(.l,T)-^(>'(.M)4  2.*K(.],2)  +  2.»K(J,3)-»H(J,4))/b. 

YT.  (T  +  nziLCTjtC'I.  Cn+7.nL(:)  +  2.»^'LC3)+f  Lr  :))/6. 

TS=YS(l4l ) 

CC  J  8  J=  J 

18  T(J)=YCJ,l4t) 

7l  =  YL(l4n 

20  TVsTv+H 

WRlTECl'^,6r:)TV.YS(,1E‘T),VtTC,JFNr)  //L(L)h::C) 

CLCnE(U:UTs9) 

CLC3E  C!:':iT=tr:) 

S7CP 

E';c 


Figure  C-2.  Program  Listing  for  1-D  Runge-Kutta  Model. 


♦  i-L  FTnTIF  TiFFFKEt.Ct  /(PFFCi' 1“ A'^ 

♦  5  MINI'IF  TIFF  STEP 
»  5  LAYFhP  Th  a 

»  IF  I,«YEPS  IF  V 

♦  UPhTM-SPuTn  LFIFt.TATiri. 

riWt*'STuM  TO,17),XKf7,lT).llK(T,ni 
AiFhAl=2.1F»l0.»* 1-7.C) 

Ai.PrtA2  =  3.»F»lo.**  1-7. 

ALFHA3=1.o7»1u.*«(-5.0} 

ALFriA<i  =  7.3F»ly.¥<(,-7.w1 
Ai,FnA5  =  6.0F»10.»*  i-7.C’» 

Sy=5.o697#ic.»*(-d.o) 

TK=219Pb. 

AlX  =  ALPhAl/t,X 

A.;X  =  ALPnA2/».X 

A3XsALPhA3/LX 

A^XsALTiiA-l/tX 

ASX=AlPHA6/LX»*2. 

ri=iAB.o?5/3./;«i)/15. 

AlYsALPKAl/LY 
AAY  =  Ai.PHA-^/LV 
A3Y=ALPnA3/CY 
A^YsALPiiA-i/i-Y 

A5Y  =  Ai,PnA5/CY*»2. 

P=3yC. 

ThFlA=(9P./2  3  77i;.)*(TF.-2)77b.) 
nPETA  =  h*9P./A3770. 
rc  3  1=1,7 
Cl  5  <^  =  1x11  4-. 

^  rFFi'(iwiT5e!FlT  t-  'F.FAThP.rAT'  ,trKM='Frt-»'ATTFC'  ,.FlA'I('5=;Cir  ;  > 
CFEh^UNlTsl6.FTLF='NSaEKP.CAT'  ,  FCKM=  •  r  OhP  AT'lE  L  '  ,  SI  A7Ut=  '  N  E  V  ) 

F90  f[pJaT{x',^°S-S  bPICGF  SUPF'CF  TtFFFKATUdFS  -  b  MIN  IS  -  5  X  Id  LAY 

l^NS'/) 

Cl  30  K=2,5lS 

9P  F£AL(S,406)i'«l,VA,TA,5F,WY 
400  FCPMAKF7.0,Fl?.3,3Fi5.t1 
IFCIMI.EO.TMI  GO  iC  1 
Go  TO  99 

'  k£xiAiX*W,'iA2X*CkY-SG*fT(l,oI  +  27  3.15)**4.)-*lAjX  +  A4X*VAW 

^WLXsuk+MXFVAJ^l-TA-T  (7,  jlj*H2A*SG»aTA-»27j.  151**4. 
1-CTC7,J)-*'2'’3.1F)**4.1 
XNfl,J)=4.*A5A*ri*(Tc2,Jj-lf l,u) J+2.*H*WSX 

XKr2,o7  =  AbX*H*fTf  3,o)-3.*Tf  2  ,uI-*2.*T  (1  .jl  1 

-.0  xk(I,Jl  =  ADX*H*fTri+l,J1-2.*Tfl,J''4TCT-1,,n'' 

Xxfo, J)=AbX*H*f2.*Tt7,.U-3.*T(fi,JJ*T (b,ul J 
<0  XN(7,J7=4.*AFX*H*(T(fi,in-l57,o))+2.*F*NLX 

Wf<VsCCSCfTHETA)*V.N/STNC  (ThFIA) 

’’C  2*=  T  =  J  ,7 

*"**u  '  -P 

WSYsAlYSwtJvix+CASY-tMYFVA)*  (Ta-T  (1, 1)  1*a2Y*SG*(  ITa  + 273.1' )**'!. 
l-(TtI, 1 )+273. 151**4.1 
Xs  1 

T^  fTVfTA  1  T  Q(j  1  XsC 

WLY  =  AlY*WNvJx+(A3Y*A4Y*VAl*CTA-T(i,l7l1-*A2Y*S0*l('lA  +  273.15l**4. 

1-(T(I, 171+273. 151**4.1 

Yn  (I,  il  =  4.*A5Y*H*CT  11,21-1(1, 11 1*2.*1;;*WSY 
YNfl,2)  =  A5Y*H*fTfI,3)-?.*l'^l,2l  +  2.*l(l,111 
CC  24  J=3,1b 

YAfi, Ol  =  AbY*M*fTfi, 0+11-2. *ia,w1+TtT,. ’-Ill 
Y^(1,161=A5Y*^*^.^.*T(1,17}- J.*YIT.  Ibl  +  nr  tbil 

25  Ya  (1  ,i7)  =  4.»  AbV*P*  Ci  (1  ,  IFl-T  (1  ,  Pi  1+2  .*li*V.T  i 
rc  Xf  T=t,7 

rC  ^5  .1=1,17 

26  TCl.Jl=I(i,o1  +  X(\(l,Jl  +  YKd,0l 
TysTW  +  n 

ThEl»  =  ThETA  +  PIHtTA  . . 

30  VnTTF(l  0  ,oOC1  X”  j’'  Cl  ,ll,T(l,2lPCl,31,Tfl,41,TCl,''),1fl,9l, 

FLPf.AlfF-^.rJi  iFo.2) 

CLCi>FtTIhIT=§) 

CLCoFCtlNTTsi^) 

51Ct> 

Fur 

Figure  C-3.  Program  Listing  for  2-D  Finite-Difference  Model. 


PM  Redacted 
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